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Abstract 

Generalized matrix Lotka-Volterra lattice equations are obtained in a systematic 
way from a "master equation" possessing a bicomplex formulation. 
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1 Introduction 

A bicomplex is an N -graded linear space (over E) M. = © S>0 A4 S together with two 
linear maps V, D : M s — > M s+1 satisfying (V - AD) 2 = for all A e R. If this holds 
as a consequence of a certain (partial differential or difference) equation, we speak of 
a bicomplex formulation for this equation 0, @ • It may be regarded as a special case 
of a parameter-dependent zero curvature condition. In this work we derive generalized 
Lotka-Volterra (LV) lattices (see ^ ||, [/J, for example) from a "master equation" 
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with a bicomplex formulation. Some consequences of the latter (conservation laws 
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Biicklund transformations 0]) are briefly discussed. We choose the bicomplex space as 
M. = Ai° ® A where A = ® s=0 A s is the exterior algebra of a 2-dimensional vector space 
with a basis £, of A 1 . It is sufficient to define bicomplex maps on since by linearity 
and action on the coefficients of monomials in £, 7] they extend to the whole of Ai (cf 
|,|). On M° we define 

5z= (2-5^)e + ir/, d^ = (S q z- z)£+ (S p z- z)r] (1) 

where i = dz/dt, S is an invertible linear operator on A4° with S — 0, and p,g G I 
The bicomplex equations for (A4, 5, d) are then identically satisfied. In the following, the 
maps 5 and d will be "dressed" , resulting in nontrivial bicomplexes. 

2 A class of generalized Lotka-Volterra equations and 
its properties 
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2.1 A bicomplex for generalized LV equations 

Let us replace d by a map D which acts as follows onM , 

Dz={QS q z-z)£+{PS p z-z) V (2) 

where Q and P are linear operators on Then (A4, 5, D) is a bicomplex iff 

QS q PS p = PS p QS q , QS q = PS q - S q - p PS p . (3) 

Choosing Ai° as the set of C m - valued functions (on an infinite lattice) z n (t), n e Z, and 
{Sz) n = Zn+i, {Qz) n = Q n z n , (Pz) n = P n z n , (|) becomes 

QnPn+q P-nQ n+p i Qn P-n Pn+q—p ■ (4) 

The ansatz 

Qn — 9 n 1 9n+q > Pn = 9n ^ 9n+p (5) 

with invertible matrices g n now solves the first equation and turns the second into 

(<?n 1 9n+q) = 9 n 1 9n+p ~ 9 n +q-p 9n+q ■ (6) 

In the following we exclude the trivial cases where p = or g = 0. If p and q are not 
coprime, let s denote their greatest common divisor, so that p = sp' and q = sq' where p', q' 
are coprime. The above equation then decomposes into independent equivalent equations 
on s sublattices. Therefore we restrict q and p to be coprime. Since {p, q, n, g n } i— > 
{— p, —q, —n,g_ n } is a symmetry of it is sufficient to consider q > 0. Generalized LV 
equations are now obtained as follows. Introducing 



V n -= g n g n +\ (7) 



equation (0) takes the form 



/t/ t/ v _ J Ki ' ' ' Vn+p-l — Vn+q-p ' " ' if > , - 

^n-"-Vn+g-lJ ~i y-1 . . . y-1 _ y-1 . . . y-1 :f < n 1°) 

I v n-l v n+p v n+q-p-l v n+q "P^-" 

In the case where 14 are scalars, we write a n instead of V n and the last two equations 
turn out to be the "extended Lotka-Volterra equations" studied in 0, §]. For q — 1, they 
reduce to <2 n = a n • ■ ■ a n _|_p_i a n _p_j_i • • ■ a n —i) and = o> n _i ■ ■ ■ fl n +p ^n+i ' ' ' ®n—pi 
respectively, which have been explored in particular in H [7|. 

Remark. In case of a finite open lattice of iV points, choose S 1 = X^i 1 ^M+i where 
are the elementary matrices with a 1 at the ith row and jth column and otherwise zeros, 
and define S 1-1 as the transpose of S. Since [S, S T ] ^ 0, (P does not yield a bicomplex. 
But after dressing, a bicomplex (yM, 5, D) for the finite open LV equation is obtained. For 



a periodic lattice choose S = J^^ 1 Ei,i+i + E 
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2.2 The bicomplex linear system and conservation laws 

The bicomplex linear system [jl] for the bicomplex introduced above is 5\ = A D% with an 
m x m matrix field x- The integrability condition 5Dx = leads to 

(Q n Xn+q Xn) Pn Xn+p Xn q-p Xn+q Xn+q-p ) (9) 

which has the form of a conservation law. As a consequence, 

+oo 

Qi = A 2^ (Qi+k(l-'P)Xi+q+k(q-p) ~ Xi+k{q-p)) (10) 
k=— oo 

where % = 0, 1, . . . , \q — p\ — 1, are conserved charges, i.e. Qi = 0. The linear system reads 
Xn - Xn+q-p = A (Q n Xn+q - Xn) , Xn = A (P n Xn+p - Xn) ■ These equations can be solved 
iteratively with a power series expansion x = J27=o Xr X (r) - Starting with x (0) = I, the 
unit matrix, we find the first set of conserved quantities 



hoo 



Q?= E (Qi+xh-p) - J ) • ( n ) 



fc=— oo 



Solving the linear system for leads to xffl ~ Xn+q-p = Qn — I and thus 



oo 



(1) -E(/-Qn-^- p) ) (12) 



with the help of which one obtains the second set of conserved charges, 

+oo oo oo 

] = (Qi+HQ-P) ~~ Qi+q+{k-j){q-p)) ~ ~ Qi+(k-j)(q-p))j (13) 

k=—oo j=l j=l 

and so forth. Expressing Q n in terms of V n via (|) and (0), conserved charges are obtained 
for the LV equations (H). 



2.3 Backlund transformations 

A "Darboux-Backlund transformation" (DBT) of a bicomplex leads to a Backlund trans- 
formation (BT) of the associated equation ||. In order to determine a (primary) auto- 
DBT of the bicomplex associated with (Q), one has to find an operator R 2 \ on Ai° which 
solves the equations [5, R21) = D 2 — Di and D 2 R21 = R21 Di where Di, D 2 are the bicom- 
plex map D depending on solutions g\ and g 2 of (0), respectively. With R 2 \ = r S p , where 
r is a function, this leads to 

r n — "Tn+q-p = Q2,n~Ql,n, ^n = P%,n — P\,n (14) 

2,n 1"n+q '"Vi Ql,n+p > p2,n ^n+p fn P\,n+p ■ (1^) 
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Using (§), (jT5|) becomes = A n = A„ +p where A n = g 2 , n r n 9i, n + P - Since are 

relatively prime, according to the Euclidean algorithm there are integers k, I such that 
kq-lp= 1. Hence A n+l = A n+1+lp = A n+kq = A n and thus 

r n = g2,lA(t)g hn+p (16) 

with an arbitrary matrix A depending on t only. ( |14"D now produces a BT: 

92,n A 9l,n+p ~ 92,n+q-p A 9l,n+q = 92,n 92,n+q ~ 9l,n 9l,n+q (17) 
(#2,tj ^<7l,n+p) = 92,n 92,n+p — 9l,n 9l,n+P ■ (18) 

In the case under consideration, the "permutability theorem" takes the form R31 + i?io = 
R32 + R20 (see @ for details). This leads to the superposition formula 

9s,n = (Ax g 2 ,n+ P - A 2 g 1>n+p ) g Q ) l+p (g^ n A 1 - g^ n Az)' 1 . (19) 

In the scalar case, the above BT can be transformed into Hirota's bilinear form [Q]. BTs 
of the scalar generalized LV equations have been derived previously in ||. 



2.4 Another class of generalized LV equations 

Instead of d, now we dress the operator 5 of the trivial bicomplex (.A/f , d, 6): 

Vz=(z- LS q - p z) £ + (i - Mz) r] (20) 

with linear operators L and M. The bicomplex conditions for (Ai,T>,d) are 

S q M - MS q = S p LS q - p - LS q , LS q ' p = MLS q ~ p - LS q - p M . (21) 

Choosing (Sz) n = z n+ i, (Lz) n = L n z n and (Mz) n = M n z n , this takes the form 

M n+g — M n = L n+P — L n , L n = M n L n — L n M n+q _ p . (22) 

Remark. L n = g n + P 9n+q an d M n = g n +p9n+p solve fl22|b) and (|22|a) becomes @, so that 
the bicomplex conditions for (A4,8, D) and (A4,T>, d) coincide. Then X> = g 5 g^ 1 and 
d = D <7 —1 , so that the two bicomplexes are gauge equivalent. ■ 

The first equation of (|2~2[) can be solved by setting 

M n = G n+P — G n , L n = G n+q — G n (23) 

with matrices G n (t). The second equation then reads 

{G n +q — G n ) = (G n+P — G n )(G n+q — G n ) — (G n+q — G n )(G n+q — G n+q ^ p ) . (24) 

Again we should restrict p, q to be coprime and different from zero. Since the map 
{p, q, n, G n } I— > {— p, —q, —n, G_„,} is a symmetry of (|24f) , it is sufficient to consider q > 0. 
Moreover, the substitution p = q — p' turns (H) with p < into the same equation where 
G n is replaced by G^ and p replaced by p'. Since > q, we thus only need to consider 
p, q > 0. Introducing W n := G n+ \ — G n we obtain 

g-l p-l 5-I g-l p-q 

- Wn+i = ( E w ^i) E - ( E ww) E w ^ • ( 25 ) 

i=0 j=0 i=0 i=o j=i-q 

The scalar case with q = 1 and p > 1 has been studied by various authors 0, |^, 0]. 
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3 Discrete time generalized Lotka-Volterra equations 

A. In order to obtain discrete time counterparts (see |7|, ||) of the generalized LV equations 
in section 2, we simply have to replace the trivial bicomplex (A4,8, d) introduced above 
by a time-discretized version: 

5z={z- TS q - p z) £ + h- 1 (Tz -z)rj, dz = {S q z - z) £ + {S p z - z)r\ (26) 

where M° is now chosen as the set of C m - valued functions z n (t) depending on two discrete 
variables n G Z, t G KL with h G R \ {0}. The operators S and T act on Ai° according 
to (Sz) n (t) = z n+ i(t) and (Tz) n (t) = z n (t + h) =: 5„(t). Next we dress d as in @ with 
(Q^)n — z n and (Pz) n = P n z n . The resulting bicomplex equations are 

Qn Pn+q Pfi Qn+p : Qn Qn h (P n Pn+q—p) ■ (27) 

The ansatz (|5]) solves the first equation and turns the second into 

9n 1 9n+q ~ 9 n X 9n+q = h (<7 n 1 g n +p — 9 n +q~p 9n+q) ■ (28) 

Again, we may restrict to coprime p, q and q > 0. Using (0), we obtain 

(I + h Vn+q-p ■ ■ ■ V n -l) V n --- Vn+q-1 = V n ■ • • (I + h V n+q ■ ■ ■ V n+p -l) (29) 

if p > q. In the scalar case with q = 1 we recover equation (3.5) in ||: 

p— 1 p— l 

«n (l + h Y\ Un-i) = «n (l + h a n+ ^j . (30) 
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If < p < q, (|D leads to 

(V n ■ ■ ■ Ki+q-p-l + hi) ~V n +q-p " ■ ■ 

= • • ■ Un.+p-i (V^+p ■ ■ • V n+q -\ + hi) (31) 

which can be mapped to (j29|) via a redefinition T i— ► TS p ~ q (so that \4(t) = \4 +p _ 9 (t + /i)) 
and a subsequent replacement (h,p, q) h- > (/i -1 , <?,£>)• Finally, if p = — r < 0, equation (|28|) 
takes the form 

9n 1 #"+9 — S'n 1 #n+g = h \{g n - r Qn) 1 ~ (dn+q Qn+q+r) ) (32) 

which in terms of V n reads 

(I + Mft^.! • • • V~ l )V n ■ ■ ■ V n+q ^ = V n ■ ■ ■ + • • • v n-r) ■ (33) 

In the scalar case with g = 1 this reduces to 

r r 



and with a change of variable v n = a„ n we recover equation (3.9) in ||. 
B. Let us now consider a dressing of 5 (while leaving d unchanged): 

Vz = (z- LTS q ~ p z)i + h- 1 {MTz- z)r] (35) 

with (Lz) n = L n z n and (Mz) n = M n z n . (A4,T>, d) is a bicomplex if and only if 

L n M n +q— p = M n L n , L n — L n+P = h (M n+g — M n ) . (36) 

Remark. Writing L n (t) = g n+p (t - h) g n+? (t) _1 and M n (t) = g n + P {t - h)g n+p {t)~ l , the 
two bicomplexes (M.,S,D) and (Ml,T>, d) are gauge equivalent. ■ 

The ansatz L n = G n — G n+q , M n = I + h {G n+P — G„) solves (|36|b) and yields the 
correct h — > limit for T>. Now (|36|a) becomes 

{G n +q — G n ) [I + h (G n+g — G n+q -p)} — [I + h (G n+P — G n )\ {G n+q — G n ) . (37) 

Again, we may assume q > 0. If p > 0, in terms of W n = G n+ \ — G n d37|) reads 

q— 1 p— q p— 1 q — 1 

( Wn +) { I+h = ( J + ^ S ^) 2 ^ • (38) 

i=0 3=1-1 3=0 i=0 

Even for g = 1 and p > 1 and in the scalar case this appears to be a new integrable time- 
discretization of the corresponding LV equation (which is a special case of (^5|)). It is 
different from the discretization (3.1) in ||. If p < 0, (^) is mapped via (p, q, h, t, G n ) \— > 
(q — p, q, —h, t + h, G^) to the above case with positive p. 



4 LV equations as reductions of a three-dimensional 
equation 

Another trivial bicomplex (Ai,5, d) is obtained by replacing (|26|) with 

5z=(z-TKA- 1 z)t + h- 1 (Tz-z)7 ] , dz = (Kz — z)£ + (Az — z)r\ (39) 

where (Tz) k>l (t) = z k> i(t + h), (Kz) k ,i = z k+1)l and (Az) k> i = z kji+x are acting on the set 
A4° of C m -valued functions z ky i(t) depending on three discrete variables k, I £ Z, t £ KL. 
The dressing Dz = (QKz -z)£+ (PAz -z)rj with Q ktl = g^j g k+u and P fc ,« = g k } gk,i+i 
then leads to a bicomplex (A4, 5, D) associated with 

9k,i - 9 k ] 9k+i,i = h {g^j g k j +1 - (}k+i,i-i 9k+i,i) ■ (40) 

In the scalar case, this equation should be equivalent to the discrete Hirota equation (cf 
@], section 4.2.5). If all the fields only depend on n := kq + Ip with fixed p,q G Z, we 
recover equation (|28|). The continuous time limit of (|40|) is 

{g k j gk+i,i)' = g k j gk,i+i - 9k+i,i ■ (41) 
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